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UV frequency metrology has been performed on the a^H - X^S+ (0,0) band of various isotopo- 
logues of CO using a frequency-quadrupled injection-seeded narrow-band pulsed Titanium:Sapphire 
laser referenced to a frequency comb laser. The band origin is determined with an accuracy of 
5 MHz {Su/u = 3 • IQ-^), while the energy differences between rotational levels in the a^II state 
are determined with an accuracy of 500 kHz. From these measurements, in combination with pre- 
viously published radiofrequency and microwave data, a new set of molecular constants is obtained 
that describes the level structure of the a^H state of ^^C^^O and '^^C^^O with improved accuracy. 
Transitions in the different isotopologues are well reproduced by scaling the molecular constants 
of ^^C^^O via the common mass-scaling rules. Only the value of the band origin could not be 
scaled, indicative of a breakdown of the Born-Oppenheimer approximation. Our analysis confirms 
the extreme sensitivity of two-photon microwave transitions between nearly-degenerate rotational 
levels of different f2-manifolds for probing a possible variation of the proton-to-electron mass ratio, 
fj, = rUp/me, on a laboratory time scale. 

PACS numbers: 33.20.-t,06.20. Jr 



I. INTRODUCTION 

The a^n state of CO is one of the most extensively 
studied triplet states of any molecule. The transitions 
connecting the a^II state to the X^S"*" ground state 
were first observed by Cameron in 1926 T. Later, the 
a?Ii state of the ^^C^^O isotopologue was studied us- 
ing radio frequency (rf) [21 [3], microwave (mw) [IHS], 
infrared [T, 'F , optical [9^ and UV spectroscopy [10] . The 
13(^16q isotopologue was studied using rf [TT] and mw [1] 
spectroscopy. 

Recently, Bethlem and Ubachs [12] identified 
metastable CO as a probe for detecting a tempo- 
ral variation of the proton-to-electron mass ratio, 
/i = mp/rrie, on a laboratory time scale. Two-photon 
microwave transitions between nearly-degenerate rota- 
tional levels in different IIsi spin-orbit manifolds were 
shown to be very sensitive to a possible variation of /i. 
As a measure of the inherent sensitivity of a transition 
to a drifting /i, the sensitivity coefficient, A'^, is defined 
via: 

^^A-,.^^ (1) 

v n 

Transitions between the J = 8, = 0, the J = 6, = 1 
and the J = A, fl = 2 levels display sensitivities ranging 
from = —300 to +200 [12] . For an overview on the 
topic of varying physical constants, we refer the reader 
to [ig and [H]. 

In this paper, we present high-precision UV measure- 
ments of the a^U - X^T,+ (0,0) band in CO. In to- 
tal 38 transitions in all six naturally occurring isotopes 
have been measured with MHz accuracy. All three ft- 
manifolds have been probed, with J up to eight. A com- 



prehensive fit of the optical data combined with previ- 
ously published rf [3] and mw [3 E] measurements was 
performed. The molecular constants found for ^^C^^O 
are mass scaled and compared with the measured transi- 
tions in other isotopologues. 

II. LEVEL STRUCTURE OF CO 

The a^H state is the first electronically excited state 
of CO, lying 6 eV above the X^T,^ ground state. CO 
in the a^II state has two unpaired electrons, leading to a 
nonzero electronic spin, S, and orbital angular momen- 
tum L. For low rotational levels, the a'^II state is best 
described in a Hund's case (a) coupling scheme, with the 
good quantum numbers A and S, the projection of L 
and S on the molecular axis, respectively. The projec- 
tions of the total angular momentum J on the molec- 
ular axis and on the space-fixed axis lead to the good 
quantum numbers and M. The basis functions are 
\nA)\v)\ST,)\jnAI), representing the electronic orbital, 
vibrational, electronic spin and rotational components of 
the wave function, respectively. For higher rotational lev- 
els, the spin decouples from the electronic angular mo- 
mentum and a Hund's case (b) coupling scheme becomes 
more appropriate. In Hund's case (b), the different fl- 
manifolds are mixed. 

The energies of the lower rotational levels of the a^H 
and X^S"*" states are shown in Fig. [l] together with 
selected transitions. The transitions are denoted by 
AJf2+i(J"), where transitions with AJ = —1,0 and 1 
are denoted by P, Q and R, respectively. As the parity 
changes in a one-photon transition, the upper lambda- 
doublet component of rotational levels in the a'^H can 
only be reached via Q-transitions, whereas the lower 
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FIG. 1. Energy level diagram of the X^E"*" {v — 0) ground 
state and the a'^II {v — 0) state of '^^C^^O. Rotational quan- 
tum numbers and total parity are listed for each level. For 
the a^Ti state, the value of the band origin, En, is subtracted 
from the energy scale. The a^Ii state has three n-manifolds, 
arising from spin-orbit coupling, and shows lambda-type dou- 
bling, as illustrated, not to scale. A number of transitions are 
indicated by the vertical arrows. 



lambda-doublet components can only be reached via P 
or R transitions. 

The separation of the electronic motion and nuclear 
motion is not exact, leading to a splitting into lambda- 
doublet states of opposite parity, as indicated, not to 
scale, in Fig. [T] The lambda doubling in the 51 = state 
is large and relatively independent of J. The lambda 
doubling in the fl = 1 and f2 = 2 manifolds is much 
smaller. In Fig. [l] the total parity, i.e., the product of 
the symmetries of the rotational and electronic parts of 
the wavefunction, is indicated by the - and + signs. The 
electronic part of the wave function of the upper (lower) 
lambda-doublet levels has / (e) symmetry. 

The spin-forbidden a^Il - system becomes 

weakly allowed due to spin-orbit mixing of singlet elec- 
tronic character into a'^II, most significantly of the A^H 
state lying 2 eV above the a^H state. As the ^^11 state 
consists of a single, f2 = 1, manifold, it only couples to 
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FIG. 2. A schematic drawing of the experimental setup. A 
pulsed beam of CO is produced by expanding CO gas into 
vacuum, using a solenoid valve. After passing through a 1 mm 
skimmer, the molecular beam is crossed at right angles with 
laser radiation at 206 nm. After being excited to the o^II 
state, the molecules fly 60 cm downstream before hitting an 
electron multiplier tube, where they are detected. The laser 
interaction region is built in a Sagnac interferometer to correct 
for Doppler shifts. 



the il = 1 levels in a'^II. Transitions to the = 0, J > 
and fl — 2 manifolds become weakly allowed by mixing of 
the different fi-manifolds. The $7 = 0, J = level is not 
mixed with the other 57-manifolds, hence, the transition 
to this level, i.e. the ^'i(l), does not obtain transition 
strength via coupling to the A^Il state. 

The a^H {v = 0) state can only decay to the ground 
state, hence, the radiative lifetimes of the different 
rotational are inversely proportional to the transition 
strengths. The lifetimes are thus strongly dependent on 
J and ri. For example the J = 2, ft ~ 2 level has a 
lifetime of 140 ms, whereas the J = 1, il = 1 level has a 
lifetime of 2.6 ms [T5] . 

Isotopes with an odd number of nucleons have a non- 
zero nuclear spin that leads to hyperfine structure. The 
two relevant odd-nucleon-number nuclei are ^'^C and "'^''O, 
with nuclear spin I — 1/2 and / — 5/2, respectively. 
Due to its zero electronic angular and spin momentum 
the hyperfine splitting in the ground-state of CO is small 
(w50 kHz) [16 . The hyperfine splittings in the a'^II state 
vary between 30 and 500 MHz for the measured transi- 
tions. 



III. EXPERIMENTAL SETUP 

The molecular beam setup used for frequency metrol- 
ogy on CO is schematically depicted in Fig. [2j A pulsed 
beam of CO is produced by expanding CO gas into vac- 
uum, using a solenoid valve (General Valve series 9). A 
backing pressure of 2 bar was used for recording transi- 
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tions at low J, while for recording transitions at higher 
J, the backing pressure was reduced to 0.5 bar. For 
recording transitions in ^■^C^^O, isotopically enriched CO 
(Linde Gas) was used. The enriched sample also con- 
tained slightly enhanced fractions of ^"^C^^O and ^■^C^''0, 
sufficient for obtaining signals. Spectra of ^^C^^O and 
12q18q TffQYQ measured using a natural CO sample. 

After passing through a 1 mm skimmer, the molecular 
beam is crossed at right angles with laser radiation tun- 
able near 206 nm. In the interaction region, a magnetic 
field of up to 200 Gauss can be applied by two coils in 
Helmholtz configuration. After being excited to the a^H 
state, the molecules fly 60 cm downstream before being 
detected by an Electron Mufiipher Tube (EMT). The re- 
sulting signal is recorded using a digital oscilloscope and 
the integrated signal is stored. The absolute detection 
efficiency of this method for metastable CO (6 eV inter- 
nal energy) is estimated to be on the order of 10~^ [T7] . 
Note that the time of flight (750 ^s) is short compared 
to the lifetime of the metastable state (>2.6 ms). An ad- 
justable slit is mounted in front of the EMT to limit the 
divergence of the beam that reaches the detector, thereby 
limiting the Doppler width of the recorded transitions. 

The spectroscopic measurements are performed with 
a narrow-band frequency-quadrupled titanium:sapphire 
(Ti:Sa) pulsed laser described in detail by Hannemann 
et al. [m HI] . A schematic drawing of the laser setup is 
shown in Fig.ls] A continuous wave (CW) Ti:Sa ring laser 
(Coherent 899) produces around 700 mW of laser power 
tunable near 824 nm. Its output is split into three parts, 
with all parts having approximately the same power. 
One part of the light is sent as a seed frequency to a 
pulsed Ti:Sa ring oscillator that is pumped at 10 Hz with 
50 mJ of pulsed 532 nm light from an injection seeded 
Nd:YAG laser (Spectra Physics LAB-170). The oscillator 
is locked to the CW seed light using a Hansch-Couillaud 
scheme. The pulsed IR light emanating from the oscil- 
lator is amplified in a bow-tie multi-pass Ti:Sa amplifier 
pumped with 300 mJ of pulsed 532 nm light from the 
same Nd:YAG laser that pumps the oscillator. After nine 
passes the laser power of the IR-beam is around 70 mJ 
in a 100 ns pulse. These pulses are then frequency dou- 
bled twice in two consecutive BBO crystals, resulting in 
pulses of 20 /iJ at 206 nm. 

In order to determine the absolute frequency of the 
CW Ti:Sa ring laser, the CW light is mixed with the light 
from an erbium-doped fiber frequency-comb laser (Menlo 
systems MComb at 250 MHz repetition frequency) that 
is locked to a global positioning system (GPS) disciplined 
Rb-clock standard. The optical interference beat signal 
is measured with a photodiode and an Agilent 53132A 
counter. The obtained beat frequency is then transferred 
via ethernet to the central computer at which the data 
is analyzed. Further details on the absolute frequency 
calibration can be found in Sec. IV Dl 

By making a small portion of the pulsed light inter- 
fere with part of the CW light, possible small differences 
between the frequencies of the CW seed laser and the 
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FIG. 3. (Color online) A schematic drawing of the laser sys- 
tem. A ring laser is used to injection seed an oscillator cavity, 
which is pumped by a NdiYAG laser at 10 Hz. The produced 
IR pulses are amplified in a multi-pass amplifier and subse- 
quently quadrupled in two consecutive BBO crystals. The 
absolute frequency of the CW seed laser is determined using 
a fiber comb laser. 01: optical isolator; GP: glass plate; OC: 
output coupler; HR: high reflective mirror; PA: Piezo Actua- 
tor; HC-lock: Hansch-Couillaud lock; SMF: single mode flber; 
AOM: acousto-optical modulator. 



central frequency of the pulsed output of the bow-tie am- 
plifier are measured and corrected for |18j . In order to 
have a good fringe visibility, and to ascertain that the 
full wavefront of the pulsed output is mapped onto the 
CW reference beam, both beams are sent through a short 
single-mode fiber. The beat pattern is detected using a 
fast photodiode in combination with an oscilloscope and 
analyzed online. 

The UV laser beam is split into two parts and sent 
through the molecular beam machine from opposite sides 
to limit the Doppler shift due to a possible imperfect 
perpendicular alignment of the laser beam. In order to 
ensure that the two beams are perfectly counterpropa- 
gating, the two laser beams are recombined after pass- 
ing through the machine, forming a Sagnac interferome- 
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FIG. 4. (Color online) Recordings of the R2{0) transition in 
i2q16q jjieasured using the laser beam propagating from the 
right hand side, upper panel, and the left hand side, lower 
panel. The frequency of the seed laser is scanned while the 
signal from the EMT is recorded. Each point represents a 
single laser pulse. For each pulse the beat frequency of both 
the fiber comb and the CW-pulse offset measurement setup 
are recorded. A single scan takes around 20 minutes. 



ter. The paths through the molecular beam machine are 
aligned such that the two beams interfere destructively 
at the exit port (a dark fringe) . The transition frequency 
is measured twice, using either the laser beam from the 
left or from the right hand side [20] . 



IV. EXPERIMENTAL RESULTS 

In Table |l] the frequencies are listed for the measured 
transitions in the a-^II - X^Y,^ (0,0) band. In view of the 
time-consuming measurement procedure, only a selection 
of transitions has been investigated. In total, 38 transi- 
tions have been recorded in all six stable isotopologues 
of CO. In i^ciSQ, 1301^0 and ^^C^^O, we have recorded 
low J transitions to each S7-manifold and each parity 
level, to be able to give a full analysis of the level struc- 
ture. Also, transitions to the J = 8, r2 = 0, J = 6, 57 = 1 
and J = 4, f2 = 2, the near degenerate levels of interest 
to the search for fi variation, were measured. Further- 
more, the i?2(0) transition, the most intense line under 
our conditions, has been measured in the ^^C^^O, ^^C^^O 
and ^'^C^^O isotopologues. Note that the observed signal 
strengths for the different transitions varies over more 
than three orders in magnitude. 

In Fig. |4] a typical recording of the i?2 (0) transition 
in ^^C^^O is shown. The upper and lower graph show 
the spectra obtained with the laser beam propagating 
through either path of the Sagnac interferometer. For 
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TABLE I. Measured transition frequencies for various iso- 
topologues of CO. The rightmost two columns are the resid- 
uals from the fit and the mass-scaling procedure discussed in 
Sec.lVll 
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FIG. 5. (Color online) Recording of the P2{7) transition 
and the Ri{7) transition in ^^C^^O. Both transitions origi- 
nate from the same ground state level and connect to one of 
two nearly degenerate levels. Therefore, the combination dif- 
ference of these two transition frequencies corresponds to the 
frequency of the transition between the two nearly-degenerate 
levels. 



each pulse of the pulsed laser, the frequency of the CW 
laser is determined using the fiber comb laser and a pos- 
sible shift between the pulsed laser and the CW laser 
is determined using the online CW-pulse offset detec- 
tion. Subsequently, these data are combined with the 
metastable CO signal from the EMT. Typically, the 
recorded scans are not perfectly linear, resulting in an un- 
even distribution of the data points along the frequency- 
axis. This has no influence on the peak determination. 
At the peak of the transition, the observed signal corre- 
sponds to typically a few thousand detected metastable 
CO molecules per laser pulse. Note that the fluctuations 
in the signal shown in Fig. |4] are due to pulse-to-pulse 
variations of the molecular beam and the UV power and 
not due to counting (Poisson) statistics. The solid lines 
in the figure show a Gaussian fit to the spectra. The tran- 
sition frequency of the i?2 (0) transition is determined by 
taking the average of the measurements taken with the 
laser beam propagating from either side. 

Transitions to the six near-degenerate levels in both 
i^C^^O and i^ciSQ were measured. For ^^C^^O transi- 
tions to four of the six near degenerate levels were ob- 
tained. In Fig. [5] a recording of the P2{7) and Ri{7) 
transitions in ^^C^^O is shown. Both transitions orig- 
inate from the J = 7 ground state level, thus the 
combination difference is equal to the frequency of the 
J = 6, ft = 1, -\ — > J = 8, ft — 0, + transition, which is 
measured to be 3256.6 MHz. 

As the transition strengths of transitions in the spin 
forbidden a^H - X^Tj^ system originate from mixing of 



FIG. 6. (Color online) A recording of the -Pi(l) transition in 
i2q16q j£ Qjjiy ^jjg coupling to the ^^11 state is considered, 
this transition has zero transition strength. 



a^H n — 1 with A^H [TS], the transition strengths of the 
different transitions are proportional to (the square of) 
the ri = 1 character of the final rotational level and the 
Honl-London factor. From this, we expect the P2{7) to be 
3.2 times more intense than the i?i(7). Experimentally, 
we find a ratio of 2.4 to 1. The deviation is explained 
by the fact that the J = 8, il — has a longer lifetime 
than the J = 6, = 1 (15.7 ms vs. 3.5 ms), and con- 
sequently, a smaller fraction of the metastable molecules 
decays back to the ground state before reaching the EMT. 
Taking the lifetime into account, we expect a ratio of 2.7 
to 1, in reasonable agreement with the experiment. 

The J = 0, f2 = is the only J = level, and is 
therefore not mixed with the 17 = 1 and fl = 2 mani- 
folds. Consequently, the i-i(l) transition, connecting the 
X^Y.+ J = 1 with the a^H fl = 0,J = 0, + parity level 
does not obtain any transition strength from coupling to 
the A^H state. Nevertheless, we were able to observe 
this transition, shown in Fig. [6] albeit with low signal 
to noise. From our measurements, we estimate that the 
Pi(l) transition is about 65 times weaker than Qi(l) and 
about 10'* times weaker than Q2{i), resulting in a life- 
time of 8(1) s [5T]. The transition strength is ascribed to 
mixing of the a^II state with a ^T,'^ state, most likely the 
ground state [22]. 

The rotational levels of isotopologues with one or two 
odd-numbered nuclei show hyperfine splitting. In Fig. [7] 
a recorded spectrum of the Q2(l) transition in -'^'^C^^O is 
shown. The transitions to the hyperfine sublcvels F = 
J - / = 1/2 and F = J + / = 3/2 of the excited state 
are clearly resolved. 
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FIG. 7. (Color online) A recording of the (52(1) transition 
in ^"^C^^O. The labels indicate the value of the total angular 
momentum F = J + / in the excited state. 



FIG. 8. The i?i(l), i?2(l) and Rail) transitions in ^^C"0 
measured in an applied magnetic field of 170 Gauss, showing 
Zeeman splitting. 



V. MEASUREMENT UNCERTAINTIES 

In this paragraph the main sources of uncertainty in 
the transition frequencies are discussed. The uncertainty 
budget is summarized in Table |Tlj 



A. Zeeman eflfect 

The ground state of CO is a ""^E state, hence, its Zee- 
man shift is small. The a^H state on the other hand has 
both electronic angular momentum and spin, and expe- 
riences a considerable Zeeman shift. In the ft = 2 state 
the effect due to the electronic angular momentum and 
due to spin will add up, whereas in the i7 = and 17 = 1 
state these effects will partly cancel. Hence, we expect 
the largest Zeeman shift to occur in the = 2 levels, in 
particular in the J = 2, = 2 level. 

In Fig. |8] the recorded spectra are shown for the 
i?2(l) and Rsi^) transitions in ^^C-'^^O in a magnetic field 
of 170 Gauss. In these measurements, the polarization 
of the 206 nm light is parallel to the applied magnetic 
field, hence, only AMj = transitions are allowed. As 
all three measured transitions originate from the A"^E+ 
J = 1 level only the Mj — ±1,0 components of the 
probed a'^II levels are observed. The Zeeman shift of 
the J = 2, 17 = 2, Mj = 1 state is 1.2 MHz/Gauss, 
while the shifts of the J ~ 2, il = 1, Mj = 1 and 
J = 2, r2 = 0, Mj — 1 are about 5 times smaller. It 
is observed that the transitions to the Mj = — 1 and 
Mj = 1 states are approximately equally strong, imply- 
ing that the polarization of the 206 nm is nearly perfectly 
linear, estimated to be better than 99%, as expected from 



light that has been quadrupled in non-linear crystals. In 
the earth magnetic field, approximately 0.5 Gauss, the 
J ^ 2, n ^ 1, Mj = -2 and Mj = 2, will be shifted 
in opposite directions by less than 1 MHz. As the po- 
larization is nearly perfect, the shift of the line center is 
estimated to be less than 1 kHz. 



B. Stark effect 

The Stark shift in the ground state of CO is neg- 
ligibly small as the dipole moment is only 0.1 Debye 
(corresponding to 0.17- 10~^cm~^/(kV/cm)) and mix- 
ing occurs between rotational levels. The Stark shift in 
the a'^H state on the other hand is larger as the dipole 
moment is 1.37 Debye and mixing occurs between the 
lambda-doublet components [55]. We estimate that the 
electric field in the excitation region is below 1 V/cm. 
This corresponds to a Stark shift of 30 kHz for the 
J = 2, il = 2, AdVl — A and less for all other levels. 



C. AC- Stark effect 

A priori, the AC-Stark shift is difficult to estimate. 
We have measured the i?2(0) transition in ^^C^^O for 
different laser powers, reducing the laser power by over an 
order of magnitude, but found no significant dependence 
of the transition frequency on laser power. Thus, we 
estimate the AC-Stark shift to be less than 100 kHz. 
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D. Uncertainty in absolute frequency 
determination 

The absolute frequency of the CW Ught is cahbrated 
by mixing this hght with the output of a frequency-comb 
laser and counting the resulting beat frequency /boat- 
The frequency /cw is then obtained by the relation 

few = n • /rep + /o + /boat; (2) 

with n the mode number of the frequency comb that is 
nearest to the frequency of the CW light, and /icp and 
fo the repetition and the carrier-envelope offset frequen- 
cies of the frequency comb, respectively, /lep is tunable 
over a small range around 250 MHz and /o is locked at 
40 MHz. We infer that the sign of /beat is positive from 
the observation that /beat increases when the CW laser 
is scanned towards higher frequency. Likewise, the sign 
of fo is positive from the observation that /beat increases 
when /o is decreased. The beat note is averaged over a 
period of 100 ms. On this time scale, the accuracy of 
the Rb-clock standard is 10"^", equivalent to 150 kHz at 
the used frequencies. This uncertainty enters separately 
in each data point taken in a frequency scan. As each 
scan is approximately 1000 data points, this uncertainty 
averages out. 

The integer mode number n is determined by measur- 
ing the transition frequency of the i?2(0) transition in 
12(-<16q ^gjjjg t]-^g frequency comb at different repetition 
frequencies and then determining at which transition fre- 
quency these three measurements coincide. This method 
gives an unambiguous transition frequency provided that 
the change in repetition frequency is much larger than 
the measurement uncertainty |24j . In our case this con- 
dition is well met as the uncertainty between two con- 
secutive measurements is on the order of 500 kHz, while 
the repetition frequency of the comb is varied by 2 MHz. 
The absolute value obtained for the -R2(0) transition of 
12(-<16q jg used to calibrate a wavelength meter (Burleigh 
WA-1500, 30 MHz precision) on a daily basis. This wave- 
length meter is then used to determine the mode number 
for the measurements of the other transitions. 



E. Doppler effect 



« I 1 ' 1 ' 1 ' 1 ' r 




Frequency - 1453233648.6 (MHz) 



FIG. 9. (Color online) The Q2(l) transition measured using 
slit widths of 10 mm and 3 mm. When the width of the slit is 
decreased the transition becomes narrower, down to 23 MHz 
for a slit width of 3 mm. The 23 MHz spectral width is 
attributed to the line width of the UV-laser. 

and 3 mm are shown. When the slit width is reduced 
below 3 mm the width of the transition remains the same 
while the signal decreases further. The minimum full 
width at half maximum (FWHM) observed is 23 MHz, 
which is attributed to the line width of the UV-radiation. 
In our measurements, we have used a slit width of 6 mm, 
as a compromise between signal intensity and line width. 

In order to eliminate the Doppler shift in the mea- 
surement, the UV laser beam is aligned in the geome- 
try of a Sagnac interferometer. The angle between the 
two counterpropagating laser beams may be estimated 
to be smaller than X/d, with A being the wavelength of 
the light in nm and d being the diameter of the laser 
beam in mm |20| . In our case this results in a maximum 
Doppler shift of 200 kHz. We have verified this by com- 
paring measurements of the i?2(0) transition of ^^C^^O 
in a pure beam of CO (velocity of 800 m/s) and a beam of 
5% CO seeded in He (longitudinal velocity of 1500 m/s). 
The second-order Doppler shift is sub-kHz. The recoil 
shift is sa 25 kHz. 



In a molecular beam experiment, the first-order 
Doppler effect is reduced by aligning the laser-beam per- 
pendicular to the molecular beam. Two residual ef- 
fects remain: (i) The finite transverse temperature of the 
molecular beam leads to a broadening of the transition, 
(ii) A possible imperfect perpendicular alignment of the 
laser beam leads to a shift of the center frequency. 

To limit the Doppler width, we have placed a variable 
slit in front of the EMT, which limits the divergence of 
the beam hitting the detector. In Fig. [9] recordings of the 
(52(1) transition of ^^C^'^^O using a slit width of 10 mm 



F. Uncertainty in peak determination 

The number of metastable molecules detected (at the 
peak of) the -R2(0) transition of ^^C'^^O, the strongest 
transition observed, is on the order of 10^ per laser pulse, 
while at the Pi(l) transition of ^^C^^O, the weakest tran- 
sition observed, it is only a few molecules per laser pulse. 
The ability to determine the line center from the mea- 
surement is limited by pulse-to-pulse variations of the 
molecular beam and the UV power. Typically, the un- 
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certainty in the peak determination is 300 kHz (corre- 
sponding to 1% of the line- width). 



G. Frequency chirp in the pulsed laser system 

The frequency of the pulsed IR light differs slightly 
from the frequency of the CW seed laser mainly due to 
two effects: (i) Mode pulling in the oscillator and (ii) Fre- 
quency chirp and shift in the amplifier. The pump laser 
induces a change of the refractive index of the Ti:Sa crys- 
tal and therefore a variation of the optical path length of 
the cavity on a time scale much shorter than the response 
time of the electronic system used for locking the cavity. 
Therefore, at the instant the pulse is produced the cav- 
ity resonances are shifted from the frequency they are 
locked to. The size of the frequency shift due to mode 
pulling depends on the pump power. The optical prop- 
erties of the Ti:Sa crystal in the bow-tie amplifier change 
depending on the population inversion, which decreases 
during the multi-step amplification of the IR pulse. The 
combined effect of both phenomena results in a shift on 
the order of a few MHz. We measure and compensate 
for these effects by making a part of the pulsed light in- 
terfere with part of the CW light. The CW-pulse offset 
can be determined within 200 kHz. However, as a result 
of wavefront distortions of the pulsed beam, it depends 
very critically on the alignment. We have measured the 
chirp at different positions in the wavefront by moving 
a small pinhole through the pulsed beam, and observed 
deviations of a few MHz. These deviations are ampli- 
fied in the harmonic generation stages. As it is unclear 
which part of the pulsed laser beam gives rise to the ob- 
served metastable CO signal, we cannot compensate for 
this effect. Consequently, we find rather large day-to-day 
deviations in the measurements; measured transition fre- 
quencies taken on the same day, with a specific alignment 
and setting of the Ti:Sa oscillator and bow-tie amplifier, 
agree within a 500 kHz. Measured transition frequencies 
taken on different days, on the other hand, may devi- 
ate by a few MHz. In order to compensate for this sys- 
tematic effect, we chose the i?2(0) transition in ^^C^^O 
as an anchor for each measurement session. With this 
anchoring, measurements of any given transition (other 
than the R2{0) transition) taken on different days agree 
within 500 kHz, showing that the consistency of the posi- 
tions of the rotational levels of a^H is sub- MHz. The rel- 
atively large systematic uncertainty due to mode pulling 
and chirp will only enter in the value of the band origin. 
The root-mean-square deviation of 14 measurements of 
the i^ciSQ i?2(0) transition is about 2 MHz. We have 
set the uncertainty of the i?2(0) transition, and hence the 
systematic uncertainty on all transitions, conservatively 
at 5 MHz. Note that Salumbides et al. have circum- 
vented this systematic offset of the absolute calibration 
by putting a small pinhole in the pulsed beam, therewith 
selecting a portion of the laser beam with smaller wave- 
front distortions. We could not apply this method here. 



Source 



Uncertainty(MHz) 



Relative: 

Zeeman effect <0.1 

Stark effect <0.03 

AC-Stark effect <0.1 

Comb absolute frequency determination <c:0.15 

Doppler effect <0.2 

Peak determination 0.3 

CW-pulse offset 0.2 
Absolute: 

CW-pulse offset 5 

TABLE II. The uncertainty budget of the measured transi- 
tions. 



as it resulted in a too large decrease in signal. 



VI. ANALYSIS 

A. Effective Hamiltonian and least square fitting 
for i^ci^o and ^^C^^O 

The effective Hamiltonian for a ^H state has been de- 
rived by several authors [101 [IS] . We have used the effec- 
tive Hamiltonian from Field et al. [TCT with the additions 
and corrections discussed by Carballo et al. [5j. The ma- 
trix elements for this effective Hamiltonian are listed in 



Table HI As discussed by Carballo et al. [5l, this Hamil- 



tonian is equivalent to the effective Hamiltonian derived 
by Brown and Merer [55], but the molecular constants 
used in these Hamiltonians have a slightly different phys- 
ical meaning, which will have consequences for the mass 
scaling discussed in Sec. |VI B} In Table [TV] relations be- 
tween the constants in the Hamiltonian of Brown and 
Merer and those in the Hamiltonian of Field et al. are 
listed for clarity. 

A least-squares fitting routine was written in Mathe- 
matica to obtain the molecular constants of the effective 
Hamiltonian. We have verified that our fitting routine 
exactly reproduces the results of Carballo et al. [5 and 
that it is consistent with PGopher |17j- For ^^C^'^O, we 
have fitted our optical data simultaneously with lambda- 
doubling transition frequencies in the rf domain mea- 
sured by Wicke et al. [3] and the rotational transition 
frequencies in the mw domain measured by Carballo et 
al. |5] and Wada and Kanamori [6]. The fitted set con- 
sists of 9 rf transitions, 28 mw transitions and 14 optical 
transitions. The different data sets were given a weight 
of one over the square of the measurement uncertainties, 
taken as 50 kHz for the mw and rf data and 1 MHz for 
the optical data. The molecular constants for the A^I]+ 
state of ^^C-'^^O are taken from Winnewisser et al. |28j . 
As discussed by Carballo et al. [5], 7 and Aj can not be 
determined simultaneously from data of a single isotopo- 
logue, therefore 7 was fixed to zero. The first column of 
Table W\ lists the different constants obtained from our 
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1) - D{x^ + 4x + 1) - A - 2Aj(i + 1) - C T C6 - 7 + a;Bo + (1 T l)(4Bi+ - 2B(J 
1) - d{x'^ + 6a; - 3) + 2C - 27 - 2B(^ + 4B+ 



En + B{x 
En + b{x 

En + 5(2: - 3) - D{x'^ - 4x + 5) + A + 2^j(i - 3) - C - 7 + - 2) 

-V2a; [B - 2D{x + 1) - - 0.67 + (1 =F "2)3+] 
-^x{x^2) [2D±B+] 

-^2{x-2) [B - 2dIx - 1) + 4j - 0.57 + B+] 



TABLE III. The matrix form of the Hamiltonian used in the comprehensive fit. x — J{J + 1). Upper sign choice refers to e 
levels, lower sign choice to / levels. 



Carballo et al. 


Brown and Merer 


En 


T + B - D - q/2 


B 


B - 2D - q/2 


D 


D 


A 


A + Ad +-f + p/2 


Aj 


Ad/2 


C 


-2/3A 


Cs 





7 


7 + P/2 


BcT 


q/2 


Bt 


p/4 + q/2 



TABLE IV. The conversion factors between the molecular 
constants used in the effective Hamiltonian of Carballo et 
al. [3 and Brown and Merer |26| . 



fit for ^^C^^O. The deviations between the observed and 
fitted transition frequencies are listed in Table |I] 

All measured transitions could be fitted to approxi- 
mately their respective uncertainties. However, the total 
root mean square (rms) of the residuals of the rf and 
mw data is significantly increased when the optical data 
is included; Carballo reported rms residuals of 27 kHz, 
while we find rms residuals of 52 kHz. The rms of the 
residuals of the fitted ^^C^^O optical transitions is equal 
to 3.3 MHz. 

The molecular constants found from the fit to the com- 
bined data agree well with the constants found from a fit 
to the rf and mw data alone, but with largely decreased 
uncertainties. The uncertainty of A, C and C^, which 
are poorly constrained by the rotational and lambda- 
doubling transitions alone, are reduced by more than a 
factor of 10. Somewhat unexpectedly, the uncertainty of 
several other constants, including -B, are also substan- 
tially decreased by the fit to the combined data. This 
can be understood from the fact that the mw data does 
not constrain B, but rather a combination of A and B. 
This also explains why the uncertainty of B as obtained 
from a fit to the rotational and lambda-doubling tran- 
sitions is ^300 kHz, whereas the rotational transitions 
have a quoted uncertainty of 5-8 kHz and are fitted with 
an rms uncertainty of 27 kHz [5]. Our optical data di- 
rectly probe A, and the more precise value of A results 
in turn in a more precise value of B. Adding the opti- 
cal data results in an uncertainty in B of 30 kHz, much 
closer to the value one would expect from the precision 



of the recorded mw transitions. Altogether, the new set 
of constants is more balanced and adequately describes 
the a?IV state. 

As seen from Table |lj the residuals of the optical tran- 
sitions probing the $7 = 1 manifold are larger than those 
to the £7 = and 17 = 2 manifolds, which is surprising 
as these transitions are the strongest transitions in the 
spectra and are measured with a higher signal-to-noise 
ratio than the other transitions. We have investigated 
whether this might be explained by perturbations arising 
from the a''^I]+ or D^/S. state. When perturbations with 
the a''^I]+ state were included, using the perturbation pa- 
rameters from Carballo et al. [S] , the residuals decreased 
marginally. A slight improvement was obtained by in- 
cluding a perturbation with the D^A state. However, as 
the perturbation parameters of the coupling between the 
a?I{ and D^A states are unknown, it is unclear if this 
improvement is genuine. 

A similar analysis has been performed for ^'^C^^O. We 
have included the hyperfine interaction in the effective 
Hamiltonian following Brown et al. |29| . Only terms that 
are diagonal in J were included, since contributions from 
off-diagonal terms are estimated to be smaller than 100 
kHz TT] . We have fitted our optical data simultaneously 
with lambda-doublet transition frequencies measured by 
Gammon et al. llj and the rotational transition frequen- 
cies measured by Saykally et al. |4j. The fitted set con- 
sisted of 19 rf transitions, 4 mw transitions and 24 optical 
transitions. The molecular constants for the X^E"*" are 
taken from Klapper et al. [TBj . The molecular constants 
resulting from the fit are listed in Table [V] The difference 
between the observed transition frequencies and the fre- 
quencies from the fit are listed in Table |l] The rms of the 
residuals of the fitted ^^C^^O optical transitions is equal 
to 0.9 MHz. 

For ^^C^^O, as for the other isotopologues, no previous 
measurements on the a^H were found in the literature, 
except for four mw-transitions in ^^C^^O f7. Hence, no 
fit has been attempted. 



B. Mass scaling 

An important motivation for this work was to validate 
the mass scaling of the energy levels of the a"^H state, 
and to confirm the sensitivity to a possible variation of 
the proton-to-electron mass ratio for a selected number 
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Molecular 












Constant 




Fitted 


Fitted 


bcaloa 


Bcaled 


E 





1453190243.1(8) 


1453340486.4(7) 


1453340489.2(9) 


1453500584(3) 


B 


-1 


50414.24(3) 


48198.28(7) 


48197.73 


45797.49 


D 


-2 


0.1919(3) 


0.1861(13) 


0.1753 


0.1582 


A 





1242751.3(10) 


1242807.6(10) 


1242806.1 


1242866.7 


Aj 


-1 


-5.732(8) 


-5.51(2) 


-5.479 


-5.206 


C 





-538.4(6) 


-536.8(3) 


-538.8 


-539.5 


Cs 





26040.4(14) 


26042.6(11) 


26044.6 


26049.4 


7 


-1 
















-2 


0.841(9) 


0.775(7) 


0.768 


0.693 




_2 




S7 81 (d\ 


o ( . y± 


ou.uo 


a 






162.2(3) 


161.9(10) 


161.9 


b 






638.0(6) 


638(2) 


638 


c 






8.3(3) 


8.5(10) 


8.5 


d 






107.3(10) 


105(3) 


105 



TABLE V. Molecular constants for the a'^Tl state of ^^C^^O and ^^C^^'O, obtained from a simultaneous fit to rf, mw and optical 
data. 7 is fixed at zero in the analysis. The one standard deviation uncertainty from the fits is listed after each value, in 
units of the last digit. The values for the band origins have an additional uncertainty of 5 MHz due to the CW-pulse offset in 
the present optical measurements. For ^'^C^^O and ^"^C^^O, the table also lists the values that were obtained by scaling the 
constants of ^^C^^O. For ^'^C'^^'O, when the constants are scaled, the hyperfine constants are fitted together with the value of 
the band origin, while the scaled constants are kept fixed, resulting in slightly difi'erent values for these constants. For ^"^C^^O 
the hyperfine parameters of ^"^C^^O are used. 



of level spittings. In the literature, the reduced mass 
of the molecule is frequently denoted by the symbol /i. 
In this paper, we will use to denote the proton-to- 
electron mass ratio, and will denote the reduced mass of 
the molecule by fired- As we will see in Sec. VIC fired is 
linearly proportional to fi which is defined for the various 
isotopologues "^C^O as 



l^red 



(3) 



The molecular constants as determined from the fits are 
effective molecular constants for the v = Q level of the 
o'^n state. In general, an effective molecular constant 
y can be expressed as 



X, 



Xe + axa(« + l/2) 



+ ax,2(f + 1/2)2 + 



(4) 



Note that for the constants A and B the second term of 
Eq. (|4| has a minus sign by convention [7 . Thus, the 
mass dependence of every constant consists of the mass 
dependence of X^ and a correction due to the vibrational 
dependence of X^^v The second column of Table [V| lists 
the dependence of the molecular constants, X^, on the 
reduced mass of the molecule, fired- The effective molec- 
ular constant, X'^^, of an isotopologue with a reduced 
mass /i^g^ then becomes 



X' — 

e,v 



^^red 
f^red 

f^'red 
f^red 

l^'red 
^red 



Xe 



ax,iiv + l/2) 



(5) 



ax,2(w + 1/2)2 + 



where X^., ax.i and ax, 2 are the constants for the iso- 
topologue with reduced mass fired and 



K. 



X 



fl dXe 
Xe dfJL 



(6) 



Note that we use fi in Eq. ^ rather than fired, see below. 

As we have only measured transitions in the v = Q 
band, the vibrational dependencies of the molecular con- 
stants cannot be extracted from our data. Hence, we have 
used the ratios between Xe, ax,i and ax, 2 determined 
by Havenith et al. [7] to scale our constants. The molec- 
ular constants for ^^C^^O and -'^'^C^^O, found by scaling 
the constants of ^'^C^'^O via the outlined procedure, are 
listed in the third and fourth column of Table |Vj respec- 
tively. For ^^C^^O, the value of the band origin and the 
hyperfine constants were determined by fitting the data 
while the other constants were fixed at the scaled val- 
ues. For ^^C^'^O, the hyperfine constants were taken to 
be identical to those of ^'^C^^O. The data for the X'^Y.+ 
state was calculated from Puzzarini et al. [30] and only 
the value of the band origin was fitted. The differences 
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FIG. 10. (Color online) In the upper panel, the measured 
values of the band origins of the a^II state of the six stable 
isotopologues of CO are plotted as a function of the reduced 
mass. The solid line shows the value of the band origin scaled 
with respect to ^^C^^O. In the lower panel the difference be- 
tween the measured and calculated values of the band origins 
are plotted. 



between the observed transition frequencies and the fre- 
quencies calculated using the scaled molecular constants 
are listed in Table |lj As is seen the correspondence is 
satisfactory. The rms of the residuals of the ^'^C^^O and 
13(^18q (j^a^a with the frequencies found by scaling the 
molecular constants is equal to 5.1 MHz and 8.3 MHz, 
respectively. 

Until now the value of the band origin. En , was treated 
as a free parameter without considering its proper mass 
scaling. The value of the band origin consists of: (i) 
A pure electronic part, that scales as (/^red)" except for 
a small correction due to the finite mass of the nuclei, 
known as the normal mass shift, or Bohr-shift, which is 
proportional to the reduced mass of the nuclei-electron 
system, (ii) A vibronic part that can be expanded in a 
power series of {v + 1/2) and contains the difference in 
zero-point energies in the a^H and X-^T,~^ states, (iii) A 



rotational part, equal to B — D, that was absorbed in 
the value of the band origin in our definition of the effec- 
tive Hamiltonian. (iv) The specific mass shift, dependent 
on the electron correlation function, (v) Nuclear-size ef- 
fects, dependent on the probability density function of 
the electrons at the nucleus. 

In the upper panel of Fig. [TO] the derived values of the 
band origin of the a^H state of the six stable isotopo- 
logues of CO are plotted as a function of the reduced 
mass. The solid line shows how the value of the band 
origin scales when effects (i-iii) , which are expected to be 
dominant, are included. The used formulas for the nor- 
mal mass shift, the vibrational and the rotational parts 
are: 



AE'n = AE„ 



AEr- 



(7) 



with 



AEnm, — E] 



no 



- 1 



(8) 




1 . . f fJred 

-- [UJenXen - WeS2^eS) —, 1 

4 \f^red 



(9) 



1 . ^ / / Mred 

8 \\^^red 



3/2 



AE^ 



f-^red 




(10) 



where AEu is the shift in the value of the band origin 
as a function of the reduced mass, /i^g^, with respect 
to a given isotopologue with reduced mass fired, band 
origin Euo and rotational constants Bq and -Do- 'niei is 
the mass of the electron. The shift was calculated with 
respect to ^^C^^O. For the Ar^E+ state the vibrational 
constants were obtained from fitting to the data from 
Coxon et al. [3T], while for the a^H state the constants 
from Havenith et al. [7| were used. In the lower panel, 
the difference between the experimental and calculated 
values of the band origins for the different isotopologues 
are plotted. The experimental and scaled values of the 
band origins deviate by a few GHz, which is more than 
expected given the precision of the data used in this anal- 
ysis. Most surprising is the difference between the val- 
ues of the band origins of ^^C^^O and ^^C^^O, which 
have nearly equal reduced masses (7.18 vs. 7.21 amu, 
respectively). This suggests a breakdown of the Born- 
Oppenheimer approximation. 
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We have re-analyzed experimental data pertaining to 
all six stable isotopologues of CO for the E^II {v = 1) 
state [32] and the C^I]+ {v — 1) state [53] and found 
a similar deviation as found for the a^II, both in size 
and in direction. This suggests that the ground state is 
probably the source of the discrepancy. Calculations us- 
ing the LEVEL program '34' were performed to estimate 
the effects of the breakdown of the Born-Oppenheimer 
approximation in the ground-state of CO following the 
approach by Coxon and Hajigeorgiou [31]. The calcu- 
lated energy difference between the ^^C^^O and ^'^C^^O 
isotopologues was approximately 50 times smaller than 
observed in our measurements and is thus insufhcient to 
explain the observed effect [35] . 

The specific mass shift was not included in our anal- 
ysis. It is however proportional to ^red, and can thus 
not explain the observed large difference between ^^C^^O 
and ^'^C^^O. Nuclear size effects could in principle cause 
a similar isotope shift as observed, but the difference in 
nuclear charge radius between ^''O and ^^O is approxi- 
mately 10 times larger than the difference between ^^O 
and ^^O, whereas the observed difference between the val- 
ues of the band origins in the isotopologues with these 
oxygen isotopes is similar in size |36| . It is therefore un- 
likely that the observed isotopic effect is due to a nuclear- 
size effect. 

C. CO(a'^n) as a target system for probing d/i/dt. 
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FIG. 11. (Color online) The frequencies of two-photon mw 
transitions between four near degeneracies as a function the 
reduced mass. The solid lines show the values obtained from 
mass scaling the molecular constants of ^^C^''0. The crosses 
indicate the calculated frequencies of the transitions in the 
six stable isotopologues, whereas the boxes show the values 
obtained from differences between measured frequencies. The 
four different transitions are listed in the legend as they ap- 
pear in the graph, from top to bottom. 



In Fig. |11[ the energies of the two-photon mw tran- 
sitions between the near degenerate levels of the differ- 
ent il-manifolds are shown as a function of ^ired- The 
solid lines show the values obtained from scaling the 
molecular constants of ^^C^^O using the mass-scaling re- 
lations discussed in Sec I VI Bl The crosses indicate the 
reduced masses of the six stable isotopologues, whereas 
the boxes show the values directly obtained from the 
measurements. The strong energy dependence on the 
reduced mass is indicative of a high sensitivity to the 
proton-to-electron mass ratio. The measured and calcu- 
lated values for the two-photon microwave transitions be- 
tween the nearly-degenerate levels are listed in Table |VlJ 
Our measurements reduce the uncertainties of the two- 
photon transitions to wl MHz, a factor 20 better than 
before. 

The mass of the proton is much larger than the masses 
of the constituent quarks, consequently, the mass of the 
proton is related to the strength of the forces between 
the quarks; Aqcdi the scale of quantum chromodynam- 
ics [33. As the same argument holds for the neutron-to- 
electron mass ratio, a possible variation of the proton- 
to-electron mass ratio is expected to be accompanied by 
a similar variation of the ncutron-to-electron mass ra- 
tio [35! • With this assumption, it follows that 



The sensitivity of a transition to a possible variation of 
the proton-to-electron mass ratio [i can now be calculated 
using the mass-scaling relations discussed before. The 
sensitivities for the two-photon microwave transitions are 
listed in the last column of Table IVII These coefficients 
have been calculated with an accuracy of 0.3-0.05%. 

It is instructive to compare the sensitivity to what one 
would expect in a pure Hund's case (a). For the tran- 
sition at 1628.3 MHz we expect a sensitivity of = 
Al2v = YimhX.ill X 1628.3 « 382 which is 12% larger 
than calculated using the model that includes the cou- 
pling between the different Jl-manifolds [35] . 

VII. CONCLUSION 

UV frequency metrology has been performed on the 
a'^H - (0,0) band of various isotopologues of CO 

using a frequency-quadrupled injection-seeded narrow- 
band pulsed Ti:Sa laser referenced to a frequency comb 
laser. 

We have fitted our optical data for ^■^C^^O together 
with the lambda-doubling transitions of Wicke et al. [5] 
and the rotational transitions of Carballo et al. [S] and 
Wada and Kanamori [6]. Adding the optical data re- 
sulted in a large decrease of the uncertainties of A, C 
and C's, and a smaller decrease in uncertainty in the other 
constants. 
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Isotopologue Transition Measured (MHz) Meas.— Calc.(MHz) K^^ 



J — 6 


Q — 
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^ ^ J — 4 


Q — 2 + 








19270.1 


3.5 


27.8 


J = b 
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, S2 = 
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1, 


> J = 4:, 


s 2 = z, — 








ibUo ( . ( 


4. / 


33. ( 


J ~ 


, i 2 = 


1 


1 \ r Q 
H ^ J = o, 










— iozo.o 


Q Q 
— O.O 


QQ /I 


T — R 

J — u 


o — 


1 

J^l 


— — r J — O, 


o — n 










/I ^ 




J = 6 


, f7 = 


1, 


+,F = 3.5- 


^ J = 4, 


= 2, 


+ ,i^ = 


6.5 


43005.8 


0.1 


12.9 


J = 6 


, n = 


1, 


-,F = 3.5 - 


^ J = 4, n 


= 2, 




6.5 


39988.0 


0.7 


12.2 


J = 6 


, n = 


1, 


+,F = 5.5- 


^ J = 8, n 


= 0, 




8.5 


22329.6 


6.0 


23.5 


J = 6 


, f7 = 


1, 


-,F = 5.5 - 


^,7 = 8, 


= 0, 


-,F = 


8.5 


4003.4 


5.0 


128 


J = 6 


, n = 


1, 


+,F = 3.5 - 


^ J = 4, n 


= 2, 


+ ,^' = 


6.5 


69062.0 


-7.1 


7.22 


J = 6 


, n = 


1, 


-,F = 3.5 - 


^ J = 4, n 


= 2, 


-,F = 


6.5 


66277.5 


-3.6 


7.60 



TABLE VI. Measured level splittings between near-degenerate levels in CO (a^H); frequencies are listed in terms of two-photon 
microwave transitions bridging the intervals. The Meas.-Calc. column lists the difference between the measured frequency and 
the frequency calculated using the constants listed in Table [Vl The sensitivities to a possible variation of the proton-to-electron 
mass ratio are listed in the last column. 



From our measurements we obtain the value of the 
band origin with an uncertainty of 5 MHz, a 30-fold 
improvement compared to the value obtained from ab- 
sorption measurements by Field et al. [TD] . We have also 
measured the value of the band origin in different isotopes 
and found an unexpected behavior of the isotope shifts, 
probably due to a breakdown of the Born-Oppenheimer 
approximation. 

Our main motivation for this study was to obtain more 
accurate values for the 2-photon transitions between near 
degenerate rotational levels in different i7-manifolds and 
validate the large sensitivity coefficients predicted for 
these transitions. The calculated values of the transitions 
agree to within a few MHz with the measured values, giv- 
ing confidence in the calculated values of K^^. 
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